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We present a measurement of the helicity of the W boson produced in top quark decays using 
ti decays in the £+jets and dilepton final states selected from a sample of 5.4 fb _1 of collisions 
recorded using the DO detector at the Fermilab Tevatron pp collider. We measure the fractions of 
longitudinal and right-handed W bosons to be fo = 0.669 ± 0.102 [±0.078 (stat.) ± 0.065 (syst.)] 
and /+ = 0.023 ±0.053 [±0.041 (stat.) ± 0.034 (syst.)], respectively. This result is consistent at the 
98% level with the standard model. A measurement with fo fixed to the value from the standard 
model yields /+ = 0.010 ± 0.037 [±0.022 (stat.) ± 0.030 (syst.). ] 

PACS numbers: 14.65.Ha, 14.70.Fm, 12.15.Ji, 12.38.Qk, 13. 38. Be, 13.88.+e 



INTRODUCTION 

The top quark is the heaviest known fundamental par- 
ticle and was discovered in 1995 [l], Q at the Tevatron 
proton-antiproton collider at Fermilab. The dominant 
top quark production mode at the Tevatron is pp — > ttX. 
Since the time of discovery, over 100 times more inte- 
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grated luminosity has been collected, providing a large 
number of ti events with which to study the properties of 
the top quark. In the standard model (SM), the branch- 
ing ratio for the top quark to decay to a W boson and a 
b quark is > 99.8%. The on-shell W boson from the top 
quark decay has three possible helicity states, and we de- 
fine the fraction of W bosons produced in these states 
as fo (longitudinal), /_ (left-handed), and f + (right- 
handed). In the SM, the top quark decays via the V — A 
charged weak current interaction, which strongly sup- 
presses right-handed W bosons and predicts fo and /_ 
at leading order in terms of the top quark mass (wit), W 
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TABLE I: Summary of the most recent W boson helicity mea- 
surements from the DO [(| and CDF Q collaborations. The 
first uncertainty is statistical and the second systematic. 

DO, f ftr 1 [6] fo = 0.425 ± 0.166 ± 0.102, 
/+ =0.119 ±0.090 ±0.053 
/+ fixed: f = 0.619 ± 0.090 ± 0.052 

fo fixed: /+ = -0.002 ± 0.047 ± 0.047 

CDF, 2.7 fb" 1 [7] f = 0.88 ± 0.11 ± 0.06, 
/+ = -0.15 ±0.07 ±0.06 
/+ fixed: f = 0.70 ± 0.07 ± 0.04 
fo fixed: /+ = -0.01 ± 0.02 ± 0.05 
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where x — Mw/m t , y — mi,/m t , and A = 1 + x A + 
y 4 — 2x 2 y 2 — 2x 2 — 2y 2 . With the present measure- 
ments of m t = 173.3 ± 1.1 GeV/c 2 3 and M w = 
80.399±0.023 GeV/c 2 [|, and taking m b to be 5 GeV/c 2 , 
the SM expected values are /o=0.698, /_=0.301, and 
/+ = 4.1 x 10~ 4 . The absolute uncertainties on the SM 
expectations, which arise from uncertainties on the par- 
ticle masses as well as contributions from higher-order 
effects, are « (0.01 - 0.02) for f and /_, and ©(lO^ 3 ) 
for /+ 1. 

In this paper, we present a measurement of the W bo- 
son helicity fractions fo and /+ and constrain the fraction 
/_ through the unitarity requirement of /_ ± / + ± fo = 1 . 
Any significant deviation from the SM expectation would 
be an indication of new physics, arising from either a de- 
viation from the expected V — A coupling of the tWb 
vertex or the presence of non-SM events in the data sam- 
ple. The most recently published results are summarized 
in Table |U 

The extraction of the W boson helicities is based on 
the measurement of the angle 6* between the opposite of 
the direction of the top quark and the direction of the 
down-type fermion (charged lepton or d, s quark) decay 
product of the W boson in the W boson rest frame. The 
dependence of the distribution of cos 9* on the W boson 
helicity fractions is given by 

w(c) cx 2(1 - c 2 )/ ± (1 - c) 2 /_ ± (1 ± c) 2 /+ (4) 

with c = cos 8*. After selection of a tt enriched sam- 
ple the four-momenta of the tt decay products in each 
event are reconstructed as described below, permitting 
the calculation of cos 9* . Once the cos 8* distribution is 
measured, the values of fo and /+ are extracted with 



a binned Poisson likelihood fit to the data. The mea- 
surement presented here is based on pp collisions at a 
center-of-mass energy y/s — 1.96 TeV corresponding to 
an integrated luminosity of 5.4 fb _1 , five times more than 
the amount used in the result in Ref. IfJ. 



DETECTOR 

The DO Run II detector [U is a multipurpose detec- 
tor which consists of three primary systems: a central 
tracking system, calorimeters, and a muon spectrometer. 
We use a standard right-handed coordinate system. The 
nominal collision point is the center of the detector with 
coordinate (0,0,0). The direction of the proton beam is 
the positive +z axis. The +x axis is horizontal, pointing 
away from the center of the Tevatron ring. The +y axis 
points vertically upwards. The polar angle, 9, is defined 
such that 9 — is the +z direction. Usually, the polar 
angle is replaced by the pseudorapidity r/ — — In tan (|). 
The azimuthal angle, </>, is defined such that = points 
along the +x axis, away from the center of the Tevatron 
ring. 

The silicon microstrip tracker (SMT) is the innermost 
part of the tracking system and has a six-barrel longi- 
tudinal structure, where each barrel consists of a set of 
four layers arranged axially around the beam pipe. A 
fifth layer of SMT sensors was installed near the beam 
pipe in 2006 @. The data set recorded before this ad- 
dition is referred to as the "Run Ha" sample, and the 
subsequent data set is referred to as the "Run lib" sam- 
ple. Radial disks are interspersed between the barrel seg- 
ments. The SMT provides a spatial resolution of approx- 
imately 10 /xm in r — <f> and 100 /im in r — z (where r is 
the radial distance in the x-y plane) and covers \r)\ < 3. 
The central fiber tracker (CFT) surrounds the SMT and 
consists of eight concentric carbon fiber barrels holding 
doublet layers of scintillating fibers (one axial and one 
small-angle stereo layer), with the outermost barrel cov- 
ering \rj\ < 1.7. The solenoid surrounds the CFT and 
provides a 2 T uniform axial magnetic field. 
The liquid-argon/uranium calorimeter system is housed 
in three cryostats, with the central calorimeter (CC) cov- 
ering \rj\ < 1.1 and two end calorimeters (EC) covering 
1.5 < 1 77 1 < 4.2. The calorimeter is made up of unit cells 
consisting of an absorber plate and a signal board; liquid 
argon, the active material of the calorimeter, fills the gap. 
The inner part of the calorimeter is the electromagnetic 
(EM) section and the outer part is the hadronic section. 
The muon system is the outermost part of the DO detec- 
tor and covers \r]\ < 2. It is primarily made of two types 
of detectors, drift tubes and scintillators, and consists of 
three layers (A,B and C). Between layer A and layer B, 
there is magnetized steel with a 1.8 T toroidal field. 
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DATA AND SIMULATION SAMPLES 

At the Tevatron, with proton and anti-proton bunches 
colliding at intervals of 396 ns, the collision rate is about 
2.5 MHz. Out of these 2.5 x 10 6 beam crossings per sec- 
ond at DO, only those that produce events which are iden- 
tified by a three-level trigger system as having properties 
matching the characteristics of physics events of interest 
are retained, at a rate of ^100 Hz [sLllol| . This analysis is 
performed using events collected with the triggers appli- 
cable for ^+jets and dilepton final states between April 
2002 and June 2009, corresponding to a total integrated 
luminosity of 5.4 fb^ 1 . Analysis of the Run Ha sample, 
which totals about 1 fb _1 , was presented in Ref. @. Here 
we describe the analysis of the Run lib data sample and 
then combine our result with the result from Ref. |6j when 
reporting our measurement from the full data sample. 

The Monte Carlo (MC) simulated samples used for 
modeling the data a re g enerated with ALPGEN in- 
terfaced to pythia [12[ for parton shower simulation, 
passed through a detailed detector simulation based on 
GEANT overlaid with data collected from a random 
subsample of beam crossings to model the effects of noise 
and multiple interactions, and reconstructed using the 
same algorithms that are used for data. For the signal 
(tt) sample, we must model the distribution of cos 6** 
corresponding to any set of values for the W boson he- 
licity fractions, a task that is complicated by the fact 
that ALPGEN can only produce linear combinations of 

V — A and V + A tWb couplings. Hence, for this analy- 
sis, we use samples that are either purely V — A or purely 

V + A, and use a reweighting procedure (described be- 
low) to form models of arbitrary helicity states. ALPGEN 
is also used for generating all F+jets processes where V 
represents the vector bosons, pythia is used for gener- 
ating diboson (WW, WZ, and ZZ) backgrounds in the 
dilepton channels. Background from multijet production 
is modeled using data. 



EVENT SELECTION 

We expect a priori that our measurement will be lim- 
ited by statistics, so our analysis strategy aims to maxi- 
mize the acceptance for tt events. The selection is done in 
two steps. In the first step, a loose initial selection using 
data quality, trigger, object identification, and kinematic 
criteria is applied to define a sample with the characteris- 
tics of tt events. Subsequently, a multivariate likelihood 
discriminant is defined to separate the tt signal from the 
background in the data. We use events in the £+jets and 
dilepton tt decay channels, which are defined below. 

In the £+jets decay tt -> W + W~bb -> Iv qq'bb, 
events contain one charged lepton (where lepton here 
refers to an electron or a muon), at least four jets with 
two of them being b quark jets, and significant miss- 



ing transverse energy $t (defined as the opposite of the 
vector sum of the transverse energies in each calorime- 
ter cell, corrected for the energy carried by identified 
muons and energy added or subtracted due to the jet 
energy calibration described below) . The event selec- 
tion requires at least four jets with transverse momen- 
tum pt > 20 GeV/c and \r]\ < 2.5 with the leading 
jet px > 40 GeV/c. At least one lepton is required 
with p T > 20 GeV/c and \r]\ < 1.1 (2.0) for electrons 
(muons). Requirements are also made on the value of 
JpT and the angle between the 1£t vector and the lep- 
ton (to reduce the contribution of events in which mis- 
measurement of the lepton energy gives rise the spurious 
$t)'- m the e+jets channel the requirement is I£t > 20 
GeV and A</)(e,^ T ) > 0.7tt - 0.045 • $t / GeV, and in 
the /i+jets channel the requirement is $t > 25 GeV and 
A<t>(H,$ T ) > 2.1 - 0.035 • $ T /GeV. In addition, for the 
/i+jets channel, the invariant mass of the selected muon 
and any other muon in the event is required to be out- 
side of the Z boson mass window (< 70 GeV/c 2 or > 100 
GeV/c 2 ). 

For the dilepton decay channel, tt — > W + W~bb — > 
ivi'v'bb, the signature is two leptons of opposite charge, 
two b quark jets, and significant J^t- The event selection 
requires at least two jets with pt > 20 GeV/c and \r]\ < 
2.5 and two leptons (electron or muon) with pt > 20 
GeV/c. The muons are required to have < 2.0, and 
the electrons are required to have |ry|<l.lorl.5<|ry|< 
2.5. 

Jets are defined using a mid-point cone algorithm 14 1 
with radius 0.5. Their energies are first calibrated to be 
equal, on average, to the sums of the energies of the par- 
ticles within the jet cone. This calibration accounts for 
the energy response of the calorimeters, the energy that 
crosses the cone boundary due to the transverse shower 
size, and the additional energy from event pileup and 
multiple pp interactions in a single beam crossing. The 
energy added to or subtracted from each jet in due to 
the above calibration is propagated to the calculation of 
$t- Subsequently, an additional correction to for the av- 
erage energy radiated by gluons outside of the jet cone 
is applied to the jet energy. Electrons are identified by 
their energy deposition and shower shape in the calorime- 
ter combined with information from the tracking system. 
Muons are identified using information from the muon 
detector and the tracking system. We require the (two) 
highest-pT lepton(s) to be isolated from other tracks and 
calorimeter energy deposits in the ^+jets (dilepton) chan- 
nel. For all channels, we require a well-reconstructed pp 
vertex (PV) with the distance in z between this vertex 
and the point of closest approach of the lepton track be- 
ing less than 1 cm. 

The main sources of background after the initial selec- 
tion in the £+jets channel are W^+jets and multijet pro- 
duction; in the dilepton channels they are Z boson and 
diboson production as well as multijet and W+jets pro- 
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duction. Events with fewer leptons than required (mul- 
tijet events, or M^+jets events in the dilepton channel) 
can enter the sample when jets are either misidentified 
as leptons or contain a lepton from semileptonic quark 
decay that passes the electron likelihood or muon isola- 
tion criterion. In all cases they are modeled using data 
with relaxed lepton identification or isolation criteria. 
The multijet contribution to the £+jets final states in 
the initially-selected sample is estimated from data fol- 
lowing the method described in Ref. [l6j |. This method 
relies on the selection of two data samples, one (the tight 
sample) with the standard lepton criteria, and the other 
(the loose sample) with relaxed isolation or identification 
criteria. The numbers of events in each sample are: 

iV loosc = N t *+ W + N M} (5) 
N tight = e,N ti+w +WV MJ (6) 

Here the coefficient eg is the efficiency for isolated 
leptons in tt or Wjjjj events to satisfy the standard 
lepton requirements, while £mj is the efficiency for a 
jet in multijet events to satisfy those requirements. 
We measure eg in Z — > tt control samples and £mj in 
multijet control samples. Inserting the measured values, 
we solve Eqs. [5] and [6] to obtain the number of multijet 
events (7V MJ ) and the number of events with isolated 
leptons (N^ +w ). In the dilepton channels we model the 
background due to jets being misidentified as isolated 
leptons using data events where both leptons have the 
same charge. This background originates from multijets 
events with two jets misidentified as leptons and from 
VF+jets events with one jet misidentified as a lepton. 

To separate the tt signal from these sources of back- 
ground, we define a multivariate likelihood and retain 
only events above a certain threshold in the value of that 
likelihood. The set of variables used in the likelihood and 
the threshold value are optimized separately for each tt 
decay channel. The first step in the optimization proce- 
dure is to identify a set of candidate variables that may 
be used in the likelihood. The set we consider is: 

• Aplanarity A., defined as 3/2 of the smallest 
eigenvalue of the normalized momentum tensor for 
the jets (in the £+jets channels) or jets and lep- 
tons (in the dilepton channels). The aplanarity A 
is a measure of the deviation from flatness of the 
event, and tt events tend to have larger values than 
background. 

• Sphericity S, defined as 3/2 of the sum of the 

two smallest eigenvalues of the normalized momen- 
tum tensor for the jets (in the ^+jets channels) or 
jets and leptons (in the dilepton channels). This 
variable is a measure of the isotropy of the energy 



flow in the event, and tt events tend to have larger 
values than background. 

• Ht, introduced in Refs. [ijj and [lH, is defined as 
the scalar sum of the jets' pt values. Jets arising 
from gluon radiation often have lower pt than jets 
in tt events, so background events tend to have 
smaller values of Ht than signal. 

• Centrality C, defined as J^- where He is the sum 
of all jet energies. The centrality C is similar to Ht 
but normalized in a way to minimize dependence 
on the top quark mass. 

• ^Tmin> defined as AR jjmin ■ ^pp , where ARj jmin 
is the distance in rj — </> space between the clos- 
est pair of jets, -Ermin is the lowest jet Et value 
in the pair, and E^f is the transverse energy of 
the leptonically-decaying W boson (in the dilepton 
channels Ej! is the magnitude of the vector sum 
of the $t and leading lepton pr). Only the four 
leading-i^ jets are considered in computing this 
variable. Jets arising from gluon radiation (as is 
the case for most of the background) tend to have 
lower values of K' Tmin . 

• mjj m i n , defined as the smallest dijet mass of pairs 
of selected jets. This variable is sensitive to gluon 
radiation and tends to be smaller for background 
than signal. 

• h, defined as the scalar sum of all the selected jet 
and lepton energies. Jets arising from gluon radia- 
tion often have lower energy than jets in tt events, 
and leptons arising from the decay of heavy flavor 
jets often have lower energy than leptons from W 
boson decay, so background events tend to have 
smaller values of h than signal. 

• Xfc: defined as the x 2 for a kinematic fit of £+jets 
final states to the tt hypothesis. Signal events tend 
to have smaller \ 2 values than background. This 
variable is not used for dilepton events, for which a 
kinematic fit is underconstrained. 

• A0(lepton, -^t), defined as the angle between the 
leading lepton and the $t- W+jets events with $t 
arising from mismeasured lepton pt tend to have 
A</>(lepton,^r) ~ or tt. 

• b jet content of the event. Due to the long life- 
time of the b quark, tracks within jets arising from 
b quarks have different properties (such as larger 
impact parameters with respect to the PV and the 
presence of secondary decay vertices) than tracks 
within light-quark or gluon jets. The consistency 
of a given jet with the hypothesis that the jet was 
produced by a b quark is quantified with a neural 
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network (NN) that considers several properties of 
the tracks contained within the jet cone 15j. In 



the £+jets channels, we take the average of the NN 
values NNj, of the two most Mike jets to form a 
variable called NNb avg , and in the dilepton chan- 
nels we take the NN& values of the two most 6-likc 
jets as separate variables NN&i (the largest NN& 
value) and NNb2 (the second- largest NN;, value). 
For top quark events, these variables tend to be 
close to one, while for events containing only light 
jets they tend to be close to zero. 

• $t or Xz- F° r the e/i and ee channels only, $t is 
considered as a variable in the likelihood discrim- 
inant. In the nfj, channel, where spurious $t can 
arise from mismeasurement of the muon momen- 
tum, we instead use Xz> the \ 2 °f a kinematic fit 
to the Z —t fifi hypothesis. 

• Dilepton mass m«. Also for the dilepton chan- 
nels only, the invariant mass of the lepton pairs is 
considered as a variable in the classical likelihood. 
The motivation is to discriminate against Z boson 
production. 

We consider all combinations of the above variables to 
select the optimal set to use for each tt decay channel. 
For a given combination of variables, the likelihood ratio 
L t is defined as 



Li 



expjEtTKDP]} 
exp{EilTH|)f]} + 



(7) 



where N val is the number of input variables used in the 
likelihood, and (Jr)^ is the ratio of the parameterized 
signal and background probability density functions. We 
consider all possible subsets of the above variables to be 
used in L t and scan across all potential selection crite- 
ria on Lt- For each Lt definition and prospective selec- 
tion criterion, we compute the following figure of merit 
(FOM): 



FOM 



y/N S + N B +<7 



(8) 



where Ns and Nb are the numbers of signal and 
background events expected to satisfy the L t selection. 

The term gb reflects the uncertainty in the background 
selection efficiency arising from any mis-modeling of the 
input variables in the MC. To assess o\b, we compare each 
variable in data and MC in background-dominated sam- 
ples. The background-dominated samples are created by 
forming a multivariate likelihood ratio (Eq. [7]) that does 
not use the variable under study, nor any variable that is 



strongly correlated with it, where the criterion is a cor- 
relation coefficient between —0.10 and 0.10. We select 
events that have low values of this likelihood, and are 
therefore unlikely to be it events, such that 95% of MC 
it events are rejected. Because the it contribution to the 
selected data sample is negligible, we can directly com- 
pare the background model to data. The impact of any 
mis-modeling on the likelihood distribution is assessed by 
taking the ratio of the observed to the expected distribu- 
tions as a function of each variable and fitting this to a 
polynomial. The result is that for each variable i we build 
a function ki that encodes the data/MC discrepancies in 
that variable. For each simulated background event, we 
reweight each likelihood according to the data/MC dif- 
ferences. For example, for a likelihood that uses n of the 
possible variables, the likelihood is given a weight 



(9) 



The quantity ob is the difference in the predicted back- 
ground yield when the unweighted and weighted L t dis- 
tributions are used for background. This uncertainty is 
propagated through the analysis as one component of the 
total uncertainty in the background yield. 



TABLE II: The set of variables chosen for use in L t for the 
e+jets and fj,+ jets channels. The numbers of background and 
tt events in the initially-selected data, as determined from a 
fit to the Lt distribution, are also presented. 





e+jets 


(U+jets 


Events passing initial selection 


1442 


1250 


Variables in best L t 


C 


C 




Ht 


Ht 










NN bavg 


NN 6avg 




Xfe 


h 










Aplanarity 




N (tt) 


592.6 ±31.8 


612.7±31.0 


N (V+jets) 


690.2 ±21.8 


579.8 ±18.6 


N (multijet) 


180.3 ±9.9 


6.5 ±4.9 



The sets of variables and L t selection criteria that max- 
imize the FOM defined in Eq. [8] for each it final state are 
shown in Tables HTland lllll Figures [T][5] show the distribu- 
tions of the variables in the best likelihood discriminant 
L t for the events passing the preselection cuts, where the 
signal and background contributions are normalized as 
described below. In addition, we use L t to determine the 
signal and background content of the initially-selected 
sample by performing a binned Poisson maximum likeli- 
hood fit to the L t distribution where the signal and to- 
tal background normalizations are free parameters. The 
W+jets contribution is determined by the fit to the L t 
distribution, while the multijet component is constrained 
to be consistent with the value determined from Eqs. [5] 
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TABLE III: The set of variables chosen for use in Lt for the dilepton channels. The number of background and tt events in the 
initially-selected data, as determined from a fit to the L t distribution, are also presented. 





e/i 


ee 


MM 


Events passing initial selection 


323 


3275 


5740 


Variables in optimized L t 




»A,tS,7Tlj j min 






K' Tmin ,$ T ,NN bl ,m u 


^ T ,NN M ,m« 


xl,NN 61 


N (tt) 


178.7 ±15.6 


74.9 ±10.7 


86.0 ±13.8 


N (background) 


144.3 ±14.5 


3200 ± 57 


5654 ± 76 



and [6] In the dilepton channels the relative contributions 
of the different background sources are fixed according 
to their expected yield, but the total background is al- 
lowed to float. The signal and background yields in the 
initially-selected sample for the £+jets channels are listed 
in Table HH and for the dilepton channels in Table IIIII 
Figures [6] and [7] show the distribution of the best like- 
lihood discriminant for each channel, where the signal 
and background contributions are normalized according 
to the values returned by the fit. Tables ITVl and W\ show 
the optimal L t cut value for each channel and the final 
number of events in data and the expected numbers of 
signal and background events after applying the Lt re- 
quirement. 



TABLE IV: Expected background and yields, and the num- 
ber of events observed, after the selection on L t in the ^±jets 
decay channels. 





e±jets 


/^±jets 


Optimized Lt 


requirement > 0.58 


> 0.29 


tt 


484.4 ±41.4 


567.2 ±47.3 


W±jets 


111. 7± 12.6 


227.7 ± 19.2 


Multijet 


58.1 ±3.9 


4.0 ±3.1 


Total 


656.2 ±43.4 


798.9 ±51.2 


Observed 


628 


803 



TABLE V: Expected background and tt yields, and the num- 
ber of events observed, after the selection on Lt in the dilepton 
decay channels. 



Source 


en 


ee 




Optimized L t requirement 


> 0.28 


> 0.934 


> 0.972 


tt 


186.6 ±0.4 


44.5 ±0.3 


43.6 ±0.3 


Z/-y* -> t+r 


N/A 


7.4 ±1.0 


19.1 ±1.3 


Z/f* — > TT 


11.2±3.7 


0.8 ±0.3 


0.35 ±0.05 


WW 


5.6 ± 1.4 


0.3 ±0.1 


0.13 ±0.05 


wz 


1.5±0.5 


0.28 ±0.04 


0.16±0.01 


zz 


1.0±0.5 


0.34 ±0.04 


0.57 ±0.04 


Misidentified jets 


15.9 ±3.1 


0.54 ±0.48 


3.7±2.5 


Total 


221.7±5.1 


54.2 ± 1.2 


67.7 ±3.9 


Observed 


193 


58 


68 



TEMPLATES 

After the final event selection, cos 6** is calculated for 
each event by using the reconstructed top quark and 
W boson four-momenta. In the £+jets decay channel, 
the four-momenta are reconstructed using a kinemetic fit 
with the constraints: (i) two jets should give the invari- 
ant mass of the W boson (80.4 GeV/c 2 ), (ii) the invariant 
mass of the lepton and neutrino should be the W boson 
mass, (iii) the mass of the reconstructed top and anti-top 
quark should be 172.5 GeV/c 2 , and (iv) the pr of the tt 
system should be opposite that of the unclustered energy 
in the event. The four highest-pr jets in each event are 
used in the fit, and among the twelve possible permuta- 
tions in the assignment of the jets to initial partons, the 
solution with the highest probability is chosen, consid- 
ering both the NN;, values of the four jets and x|. This 
procedure selects the correct jet assignment in 59% of MC 
tt events. With the jet assigned, the complete kinemat- 
ics of the tt decay products (i.e., including the neutrino) 
are determined, allowing us to boost to the rest frames of 
each W boson in the event. We compute cos 8* for the W 
boson that decays leptonically. The hadronic W boson 
decay from the other top quark in the event also contains 
information about the helicity of that W boson, but since 
we do not distinguish between jets formed from up-type 
and down-type quarks, we can not identify the down-type 
fcrmion to calculate cos 6**. We therefore calculate only 
|cos0*|, which is identical for both jets in the rest frame 
of the hadronically decaying W boson. Left-handed and 
right-handed W bosons have identical |cos^*| distribu- 
tions, but we can distinguish either of those states from 
longitudinal W bosons, thereby improving the precision 
of the measurement. 

In the dilepton decay channel, the presence of two neu- 
trinos prevents a constrained kinematic fit, but with the 
assumption that the top quark mass is 172.5 GeV/c 2 , an 
algebraic solution for the neutrino momenta can be ob- 
tained (up to a two-fold ambiguity in pairing the jets and 
lcptons, and a four- fold solution ambiguity). To account 
for the lepton and jet energy resolutions, the procedure 
described above is repeated 500 times with the energies 
fluctuated according to their uncertainties, and the aver- 
age of all the solutions is used as the value of the cos 9* for 
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FIG. 1: (Color online) Comparison of data and MC of the variables for preselected events, chosen for the best likelihood 
discriminant L t in the e+jets channel: (a) ^4, (b) C, (c) Ht, (d) xli ( e ) m jjmin, (f) K^mim an d (g) NNi, avg . The uncertainties 
on the data points are statistical only. 
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each top quark. 

As mentioned above, the extraction of both /o and 
/+ requires comparing the data with the MC models in 
which both of these values are varied. Since ALPGEN can 
only produce linear combinations of V — A and V + A 
tWb couplings, it is unable to produce non-SM / values, 
and can produce / + values only in the range [0, 0.30]. We 
therefore start with ALPGEN V — A and V + A samples, 
and divide the samples in bins of parton-level cos 6*. For 
each bin, we note the efficiency for the event to satisfy 
the event selection and the distribution of reconstructed 
cos 9* values. With this information we determine the 



expected distribution of reconstructed cos 9* values for 
any assumed W helicity fractions, and in particular we 
choose to derive the distributions expected for purely 
left-handed, longitudinal, or right-handed W boson, as 
shown in Fig. [5] The deficit of entries near cos6>* = — 1 
relative to the expectation from Eq.@]is due to the px re- 
quirement imposed when selecting leptons. We verify the 
reweighting procedure by comparing the generated V±A 
ALPGEN samples with the combination of reweighted dis- 
tributions expected for V ± A couplings, and find that 
these distributions agree within the MC statistics. The 
templates for background samples are obtained directly 
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FIG. 3: (Color online) Comparison of data and MC of the variables for preselected events, chosen for the best likelihood 
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FIG. 4: (Color online) Comparison of data and MC of the variables for preselected events, chosen for the best likelihood 
discriminant L t in the ee channel: (a) A, (b) S, (c) mjj m i n , (d) |?t, (e) NNf,i, and (f) m«. The uncertainties on the data 
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from the relevant MC or data background samples, and 
are shown in Fig. [SJ 

MODEL-INDEPENDENT W HELICITY FIT 

The W boson helicity fractions are extracted by com- 
puting a binned Poisson likelihood L(/o,/+) with the 
distribution of cos#* in the data to be consistent with 
the sum of signal and background templates. The like- 
lihood is a function of the W boson helicity fractions fo 



and /+, defined as 

Hf ,f+) = n n e-("M;-n>, y ) 2 /2< y x 

i=l j=l 

iV b i„ M 

II p ( d ik-,nik) (10) 
fe=l 

where P{dik' 1 nik) is the Poisson probability for observing 
dik events given a mean expectation value riik, -/Vchan is 
the number of channels in the fit (a maximum of five 
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in this analysis: e+jets, /i+jets, e/i, ee, and fJ.fJ.), A^kg,i 
is the number of background sources in the i th channel, 
-^bins,i is the number of bins in the cos 9* distribution 
for any given channel (plus the number of bins in the 
| cos 9* | distribution for hadronic W boson decays in the 
^+jets channels), n^y is the nominal number of cos#* 
measurements from the j th background contributing to 
the i th channel, <?h,ij is the uncertainty on Hb,ij, n>b,ij is 
the fitted number of events for this background, dik is 
the number of data events in the k th bin of cos 9* for the 
i th channel, and is the predicted sum of signal and 



background events in that bin. The can be expressed 
as 

_ eahPa.ik + e+f + p + . ik + e_(l - / - f+)p-,ik 
f-s- + / £o + /+£+ 

^Vbkg 

+ nb,ijPb,ijk (11) 

3=1 

where n Sj , represents the number of cos 9* measurements 
from signal events in a given channel, the p represent the 
probabilities for an event from some source to appear in 
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FIG. 6: (Color online) Best L t variable for the (a) /x+jets 
and (b) e+jets channels. The normalization of the signal and 
background models is determined by the Poisson maximum 
likelihood fit to the Lt distribution. The arrows mark the 
required Lt values for events in each channel. 

bin k for channel i (as determined from the templates), 
and the subscripts 0, +, — refer to the templates for 
tt events in which the W bosons have zero, negative, 
or positive helicity, and the subscript b, i refers to the 
templates for the i th background source. The efficiency 
for a event to satisfy the selection criteria depends upon 
the helicity states of the two W bosons in the event; the e 
are therefore necessary to translate the fractions of events 
with different helicity states in the selected sample to the 
fractions that were produced. The quantity E\ is defined 
as 

£A =]T/v£AA' (12) 

A' 

where £aa' is the relative efficiency for events with W 
bosons in the A and A' helicity states to satisfy the selec- 
tion criteria. The values of e\x' for each tt decay channel 
are given in Table I VII While performing the fit, both 
/o and /+ are allowed to float freely, and the measured 
W helicity fractions correspond to those leading to the 
highest likelihood value. 

We check the performance of the fit using simulated 
ensembles of events, with all values of /o and /+ from 



E (c) D0, L=4.3fb 1 -^Data 
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FIG. 7: (Color online) Best Lt variable for the (a) e/x, (b) ee 
and (c) [i{x decay channels. The normalization of the signal 
and background models is determined by the Poisson maxi- 
mum likelihood fit to the L t distribution. The arrows mark 
the required Lt values for events in each channel. 



through 1 as inputs in increments of 0.1, with the sum of 
/o and /+ not exceeding unity. We simulate input data 
distributions for the various values by combining the pure 
left-handed, longitudinal, and right-handed templates in 
the assumed proportions. In these ensembles, we draw 
a random subset of the simulated events, with the num- 
ber of events chosen in each channel fixed to the number 
observed in data. Within the constant total number of 
events, the numbers of signal and background events are 
fluctuated binomially around the expected values. Each 
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FIG. 8: Distribution of cos 9* in tt MC samples that were 
reweighted to derive the distributions for purely left-handed, 
longitudinal, or right-handed W bosons. The distribution for 
leptonically- and hadronically-decaying W bosons in ^+jets 
events are shown in (a) and (b), respectively, and the distri- 
bution for dilepton events is shown in (c). For hadronically 
decaying W bosons the cos 0* distribution for left- and right- 
handed W bosons are identical. All of the distributions are 
normalized to unity. 
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FIG. 9: (Color online) Distribution of cos 6* in background 
samples. The distribution for leptonically- and hadronically- 
decaying W bosons in £+jets events are shown in (a) and (b) , 
respectively, and the distribution for dilepton events is shown 
in (c). All of the distributions are normalized to the expected 
yield for each source of background. 



of these sets of simulated events is passed through the 
maximum likelihood fit using the standard cos 8* tem- 
plates. We find that the average fit output value is close 
to the input value across the entire range of possible val- 
ues for the helicity fractions, with the small differences 
between the input and output values being consistent 
with statistical fluctuations in the ensembles. As an 
example, the set of /o and /+ values obtained when it 
events are drawn in the proportions expected in the SM 
is shown in Fig. [TO] 



SYSTEMATIC UNCERTAINTIES 

Systematic uncertainties are evaluated using simulated 
event ensembles in which both changes in the background 
yield and changes in the shape of the cos 8* templates in 
signal and background are considered. The simulated 
samples from which the events are drawn can be either 
the nominal samples or samples in which the systematic 
effect under study has been shifted away from the nom- 
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TABLE VI: Efficiencies of different W boson helicity config- 
urations in tt events to pass the selection criteria, relative to 
the efficiencies for a mixture of V — A and V + A events. 
The indices — ,0 and + correspond to the helicity states of 
the two W bosons, and their order is leptonic W, hadronic 
W for the £+jets channel, and arbitrary for dilepton channels 
(where there is no distinction between the two W bosons in 
the event). Small differences in values in the dilepton chan- 
nels under interchange of the indices are from variations in 
MC statistics. 





e+jets /^+jets 


e[i 


ee 




£ 


0.76 


0.73 


0.67 


0.68 


0.68 


£-0 


0.87 


0.83 


0.84 


0.86 


0.85 


£-+ 


0.76 


0.73 


0.88 


0.89 


0.89 


£0- 


0.94 


0.95 


0.85 


0.86 


0.87 


£00 


1.08 


1.09 


1.06 


1.05 


1.05 


£0+ 


0.94 


0.95 


1.10 


1.05 


1.05 


£+- 


0.92 


0.96 


0.89 


0.88 


0.91 


£+0 


1.06 


1.11 


1.12 


1.03 


1.07 


£++ 


0.92 


0.96 


1.15 


0.99 


1.03 
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FIG. 10: Fit values for f and /+ obtained with 1000 MC 
simulations of the W boson helicity measurement. The SM 
helicity fractions, marked by the star, were taken as input to 
the simulations. The triangle corresponds to the physically 
allowed region where fo + f+ < 1. 



inal value. In general, the systematic uncertainties as- 
signed to /o and /+ are determined by taking an average 
of the absolute values of the differences in the average 
fit output values between the nominal and shifted V — A 
and V + A samples. 

The jet energy scale, jet energy resolution, and jet iden- 
tification efficiency each have relatively small uncertain- 
ties that are difficult to observe above fluctuations in the 
MC samples. To make the effects more visible, we vary 
these quantities by ±5 standard deviations, and then di- 



vide the resulting differences in the average fit output 
by 5. The top quark mass uncertainty corresponds to 
shifting m t by 1.4 GeV/c 2 , which is the sum in quadra- 
ture of the uncertainty on the world average m t (1.1 
GeV/c 2 ) and the difference between the world average 
value (173.3 GeV/c 2 ) and the value assumed in the anal- 
ysis (172.5 GeV/c 2 ). We evaluate the contribution of 
template statistics to the uncertainty by repeating the fit 
to the data 1000 times, fluctuating the signal and back- 
ground distributions according to their statistics in each 
fit. The uncertainties due to the modeling of tt events 
are separated into several categories and evaluated us- 
ing special-purpose MC samples. The uncertainty in the 
model of gluon radiation is assessed using PYTHIA MC 
samples in which the amount of gluon radiation is shifted 
upwards and downwards; the impact of NLO effects is as- 
sessed by comparing the default leading-order ALPGEN 
generator with the NLO generator MC@NLO [l!|; the 
uncertainty in the hadronic showering model is assessed 
by comparing ALPGEN events showered with PYTHIA and 
with HERWIG [13] ; and lastly, the impact of color recon- 
nection effects is assessed by comparing PYTHIA samples 
where the underlying event model does and does not in- 
clude color reconnection. The uncertainty due to data 
and MC differences in the background cos 0* distribution 
is derived by taking the ratio of the data and the MC dis- 
tribution for a background-enriched sample (defined by 
requiring that events have low values of L t ) and then us- 
ing that ratio to re- weight the distribution of background 
MC events that satisfy the standard selection. The un- 
certainty in the heavy flavor content of the background is 
estimated by varying the fraction of background events 
with heavy flavor jets by ±20%. Uncertainties due to 
the fragmentation of b jets are evaluated by comparing 
the default fragmentation model, the Bowler scheme 21 1 
tuned to data collected at the CERN LEP collider, with 
an alternate model tuned to data collected by the SLD 
collaboration [22| . Uncertainties in the parton distribu- 
tion functions (PDFs) are estimated using the set of 2 x 20 
errors provided for the CTEQ6M ^23] PDF. The analysis 
consistency uncertainty reflects the typical difference be- 
tween the input helicity fractions and the average output 
values observed in fits to simulated event ensembles. Fi- 
nally, we include an uncertainty corresponding to muon 
triggers and identification, as control samples indicate 
some substantial data/MC discrepancies for the loose se- 
lection we use. All the systematic uncertainties are sum- 
marized in Table IVITl 
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TABLE VII: Summary of the absolute systematic uncertain- 
ties on f + and fo- 



Source 


Uncertainty (/ + 


) Uncertainty (fo) 


Jet energy scale 


0.007 


0.009 


Jet energy resolution 


0.004 


0.009 


Jet ID 


0.004 


0.004 


Top quark mass 


0.011 


0.009 


Template statistics 


0.012 


0.023 


ti model 


0.022 


0.033 


Background model 


0.006 


0.017 


Heavy flavor fraction 


0.011 


0.026 


b fragmentation 


0.000 


0.001 


PDF 


0.000 


0.000 


Analysis consistency 


0.004 


0.006 


Muon ID 


0.003 


0.021 


Muon trigger 


0.004 


0.020 


Total 


0.032 


0.060 



RESULT 

Applying the model independent fit to the Run lib 
data, we find 

fa = 0.739 ± 0.091 (stat.) ± 0.060 (syst.) (13) 
/+ = -0.002 ±0.045 (stat.) ± 0.032 (syst.). 

The comparison between the best-fit model and the 
data is shown in Fig. QTJ and the 68% and 95% C.L. 
contours in the (/+,/o) plane are shown in Fig. IT2Ta'). 
To account for systematic uncertainties, we perform a 
MC smearing of the L distribution, where the width of 
the smearing in fo and /+ is given by the systematic 
uncertainty on each helicity fraction, and the correlation 
coefficient of —0.83 between them is taken into account. 

To assess the consistency of the result with the SM, 
we note that the change in — lnL(/o,/+) (Eq. ITU]) be- 
tween the best fit and the SM points is 0.24 considering 
only statistical uncertainties and 0.16 when systematic 
uncertainties are included. The probability of observing 
a greater deviation from the SM due to fluctuations in 
the data is 78% when only the statistical uncertainty is 
considered and 85% when both statistical and systematic 
uncertainties are considered. 

We have also split the data sample in various ways to 
check the internal consistency of the measurement. Using 
^±jets events only, we find 

fo = 0.767 ±0.117 (stat.), (14) 
/+ = 0.018 ±0.061 (stat.); 

and when using only dilepton events we find 

fa = 0.677 ±0.144 (stat.), (15) 
/+ = -0.013 ±0.065 (stat.). 

We also divide the sample into events with only 
electrons (e±jets and ee) and events with only muons 



(/x±jets and /x/x). The results for electrons only are 



fa = 0.816 ±0.142 (stat.), (16) 
/+ = -0.063 ±0.066 (stat.), 

and for muons only are 

fo = 0.618 ±0.150 (stat.), (17) 
/+ = 0.130 ±0.081 (stat.). 



Finally, we perform fits in which one of the two helicity 
fractions is fixed to its SM value. Constraining fa, we find 

/+ = 0.014 ± 0.025 ± (stat.) ± 0.028(syst.), (18) 

We also constrain / + and measure fo, finding 

fa = 0.735 ± 0.051 (stat.) ± 0.051(syst.). (19) 



COMBINATION WITH OUR PREVIOUS 
MEASUREMENT 

To combine this result with the previous measurement 
from Ref. Q , we repeat the maximum likelihood fit with 
the earlier and current data samples and their respective 
MC models, treating them as separate channels in the 
fit. This is equivalent to multiplying the two-dimensional 
likelihood distributions in fa and /+ corresponding to the 
two data sets. We determine the systematic uncertainty 
on the combined result by treating most uncertainties as 
correlated (the exception is template statistics) and prop- 
agating the uncertainties to the combined result. The 
results are presented in Table IVIIII 



TABLE VIII: Summary of the combined systematic uncer- 
tainties on /+ and fo for Run Ha and Run lib. 



Source 


Uncertainty (/+) Uncertainty (fo) 


Jet energy scale 


0.009 


0.010 


Jet energy resolution 


0.004 


0.008 


Jet ID 


0.005 


0.007 


Top mass 


0.012 


0.009 


Template statistics 


0.011 


0.021 


ti model 


0.024 


0.039 


Background model 


0.008 


0.023 


Heavy flavor fraction 


0.010 


0.022 


6 fragmentation 


0.002 


0.004 


PDF 


0.000 


0.001 


Analysis consistency 


0.004 


0.006 


Muon ID 


0.002 


0.017 


Muon trigger 


0.003 


0.024 


Total 


0.034 


0.065 



The combined result for the entire 5.4 fb 1 sample is 

fo = 0.669 ± 0.078 (stat.) ±0.065 (syst.), (20) 
/+ = 0.023 ±0.041 (stat.) ±0.034 (syst.). 
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FIG. 11: (Color online) Comparison of the cos#* distribution 
in Run lib data and the global best-fit model (solid line) 
and the SM (dashed line) for (a) leptonic W boson decays in 
^+jets events, (b) hadronic W boson decays in £+jets events, 
and (c) dilepton events. 



The combined likelihood distribution is presented in 
Figs. HUJb). The probability of observing a greater de- 
viation from the SM due to fluctuations in the data is 
83% when only statistical uncertainties are considered 
and 98% when systematic uncertainties are included. 
Constraining fo to the SM value, we find 

/+ = 0.010 ± 0.022 (stat.) ± 0.030 (syst.) (21) 

and constraining /+ to the SM value gives 

fo = 0.708 ± 0.044 (stat.) ± 0.048 (syst.). (22) 

CONCLUSION 

We have measured the helicity of W bosons arising 
from top quark decay in tt events using both the £+jets 



and dilepton decay channels and find 

fo = 0.669 ± 0.102 (23) 
[±0.078 (stat.) ± 0.065 (syst.)], 

/+ = 0.023 ±0.053 

[±0.041 (stat.) ± 0.034 (syst.)]. 

in a model-independent fit. The consistency of this mea- 
surement with the SM values / = 0.698, /+ = 3.6 x 10~ 4 
is 98%. Therefore, we report no evidence for new physics 
at the tWb decay vertex. 
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